Electron pairing in one-dimensional quasicrystals 
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Electron pairing in one-dimensional binary Hubbard chains is studied for different values of the 
band-filling using the Density Matrix Renormalization Group method. The systems consist of linear 
arrays of sites with two types of on-site correlations defined by two potentials: Ua being attractive 
(< 0), and Ub repulsive (> 0). The atomic levels of the system are modulated with periodic and 
quasiperiodic ordering, in the latter case following the Fibonacci sequence. We analyze the effect of 
such modulations and calculate the electron pairing phase diagram as a function of the band-filling. 
It is observed that there is a critical value of the band-filling where the behavior of the periodic and 
the Fibonacci binary Hubbard chains is reversed. 

PACS numbers: 71.10.Fd, 71.10.Li, 71.23.Ft 



An impression that low-dimensional systems are an 
ideal and purely serve as a mathematical example is 
vanishing from stage. The several nanoscopic low- 
dimensional devices that have been obtained in the labo- 
ratories such as carbon nanotubes, semiconducting quan- 
tum wires and quantum dots mean a breakthrough in the 
engineering of materials with novel physical properties; 
with desired properties such as superconductivity, and 
with promising applications, for example, in the area of 
spintronics. Such systems are inhomogeneous in their 
structure and intrinsically strongly correlated due to the 
reduction of the available phase space. To this collection 
of materials we can add the one-dimensional (ID) qua- 
sicrystals (QCs). One of the best known ID quasicrys- 
tals is based on the Fibonacci sequence which started 
drawing interest after the papers by Kohmoto et alji and 
Ostlund et alA. The spectral properties of the Fibonacci 
chain are exotic; the single-particle eigenstates are nei- 
ther extended nor localized but critical and the spec- 
trum corresponds to that of a Cantor setji~— . A Fibonacci 
sequence consists of two elements A and B and the en- 
tire sequence is generated by successive application of the 
substitution rule. The first few generations are Gq = B, 
d = A, G 2 = AB, G 3 = ABA, G 4 = ABAAB, 
Gi — Gi-iGi-2 for i > 2, where the letter Gi indicates 
the ith generation. In a Fibonacci chain, the elements 
A and B from the Fibonacci sequence may denote two 
different atoms (diagonal model) or two different bonds 
separating identical atoms (off-diagonal model). In this 
work, we will study the diagonal model, where the site 
energy takes two values ea and eb associated to atoms 
A and B, respectively. In the corresponding diagonal 
model the number of sites with energy ea is Na(ti) and 
the number of sites with energy eb is Nb(ji). The total 
number of sites in a generation n is represented by N(n), 
N(0) = N(l) = 1. These numbers are related by 

N(n) = N(n-l)+N(n-2), 
N A (n) = N(n-l), 
N B (n) = N(n-2). 



In the quasiperiodic limit (n — > oo) the ratio 
NA{n)/NB(n) converges towards the golden mean a = 
(V5 + l)/2. Motivated by the work of Alexandrov et al£ 
on the s— wave electron pairing in a one-dimensional bi- 
nary Hubbard system with two electrons and a total mo- 
mentum K = 0, we investigate the electron pairing for 
different values of the band-filling in a Fibonacci lattice. 
In their paper, Alexandrov et al. considered the possibil- 
ity of a binary system due to the fact that the C uO chains 
play a key role in the high-temperature superconductiv- 
ity. Thus, they considered a Hubbard model with repul- 
sion on copper and attraction on oxygen atoms. And 
for simplicity, the energy of the atomic levels was set the 
same for all sites. We further want to study the interplay 
of the binary Hubbard systems under the influence of a 
modulation of the atomic energy levels, which we model 
by setting their values along the chain distributed accord- 
ing to the Fibonacci sequence. The results are compared 
to those of the homogeneous and diatomic chains. The 
binary Hubbard system proposed by Alexandrov with 
only two electrons has been solved by means of the Bethe 
Ansatz method. In our case, we use the Density Matrix 
Renormalization Group method (DMRG) 6 ~ 9 which is an 
efficient method for investigating low-energy properties 
of many-body and strongly correlated systems such as 
those briefly described above. 

A one-dimensional binary Hubbard chain consists of 
non-equivalent sites distinguished through two on-site 
potentials: Ua being attractive (< 0), and Ub repulsive 
(> 0) arranged in alternating order. The one-dimensional 
Hubbard model is one of the few examples of an ex- 
actly solvable model using Bethe Ansatz^, where the 
./V-particle wave function is constructed using plane- wave 
exponents with coefficients obtained from a two-particle 
S-matrix. The binary Hubbard system in Ref. [H was 
formed by a collection of unit cells consisting of two sites: 
A and B (see Figure la). For the zth unit cell the cre- 
ation (annihilation) operators with spin a (=1~, i) are a\ a 
{a.i,a) and b\ a (&i i<7 ) at sites A and B, respectively. The 



following Hamiltonian was then proposed: 

H = -t^ {4,<7 b i-l,° + b \,a a i+l,<y + a \,<jK° + b l<T a i,<r} 
i,cr 

+E{^<t<+^<4}, 

i 

with ea — £b =0 and two electrons in the chain, nf a 
and n\ a are the the electron number operators on sites 
A and B, respectively, with Ua the attractive potential 
on site A and Ub the repulsive potential on site B. From 
the Bethe Ansatz solution it was then concluded that, 
given a value Ub for the repulsive interaction, there is a 
critical value of Ua above which pairing of the electrons 
was always taking place in the system. The value of the 
attractive interaction is given by: 



\Ua\> 



2U b 
U b + ' 



(1) 



In Ref. it was also left open what would happen for 
systems with more than two electrons. One possibility 
mentioned there was that doping the systems would ren- 
der in less paired-states, which would agree with the fact 
that for high-temperature superconductors there is a sat- 
uration in their critical temperature when doping these 
materials. 

We want now to extend the investigation of electron 
pairing to binary Hubbard systems with more than two 
particles and to lattices with different topology like the 
Fibonacci one. We propose the following Hamiltonian: 



H = ^2 Sim 



■ff.c.+y^ Uinifmi, (2) 



where c] CT (ci. CT ) is the creation (annihilation) operator 

with spin a (=t)4-) at site i and 7ij CT = c\ a a^ is the 
electron number operator. Ei is the atomic energy level 
and t = 1 is the nearest neighbor hopping matrix, which 
we choose to set the energy scale. Ui is the on-site in- 
teraction and will take the value Ua for i odd and U b 
for i even. The Hamiltonian in Eq. @ incorporates the 
different systems we want to study. We want to compare 
three different atomic level scenarios: i) The homoge- 
neous case, which corresponds to the system studied in 
Ref. [U, and in which = e, i.e., the atomic level is set 
the same for all sites and we set to zero, ii) The diatomic 
case, where there are two different values intercalated, Ea 
and eb { £a < £b)- In this case, half of the sites have 
an atomic level of ea and half of them have the value 
Eb- Hi) The last system is the Fibonacci case, where 
we assign each site an atomic energy value following the 
Fibonacci sequence starting with ea and thus we will 
have more sites with ea than with eb- The exact num- 
ber depends on the generation of the Fibonacci set to be 
used. Figures la, lb, and lc show the atomic level struc- 
tures just described. It is expected that the quantum 
confinement originated in the diatomic and Fibonacci 
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FIG. 1: Structure of one-dimensional binary Hubbard Sys- 
tems, a) The homogeneous binary Hubbard chain, with 
ea — £b — £- b) The diatomic binary Hubbard chain, with 
ea eb- c) The Fibonacci binary Hubbard chain, with 
ea 7^ Eb- Up to the 5th Fibonacci generation is shown in 
this diagram and ea < eb for b) and c). 



chains due to the fact that ea < Eb will, in general, 
enhance localization and therefore favor the formation of 
local electron pairs. On the other hand, the total effect 
of the different local potentials and confinement might 
be cumbersome. Furthermore, at half-filling Umklapp 
scattering plays a key role in the electronic properties 
of ID systems. As already mentioned, the interesting 
materials and systems have an inhomogeneous structure, 
which on the atomic level induces potentials that modify 
their properties. In order to investigate such effects we 
need to consider the strong correlations in many-body 
quantum systems. Lacking an analytical solution, we 
make our attempts using the density matrix renormal- 
ization group method (DMRG). This method has its ori- 
gin in the numerical renormalization group formulated by 
Wilsonii and allows for accurate calculations of ground 
state properties in low-dimensional quantum lattice sys- 
tems for which the basis of the Hilbert space grows ex- 
ponentially with the number of particles and cannot be 
handled using exact diagonalization. The DMRG is a 
variational, real-space method that selects in a system- 
atic way a sector of the Hilbert space that best repre- 
sents the ground state of a system, which is done by se- 
lecting only the m most probable states calculated from 
the density matrix of the system. The numerical error 
caused by truncation of the original basis can be calcu- 
lated directly as the total weight of the states that were 
discarded. For our systems, with L = 144 sites, we kept 
m = 256 density-matrix states, resulting in a maximum 
truncation error of the order of 10 -6 . The length of the 
systems corresponds to the 12th Fibonacci generation, 
being eb the zeroth generation and ea the first one. The 
systems above described were investigated for different 
values of the band-filling under open boundary condi- 
tions, which favors convergence in the DMRG method. 
Finite-size effects such as particle density oscillations and 
charge accumulation close to the edges of the systems are 
present specially in the homogeneous system away from 
half-filling. The diatomic and Fibonacci chains behave 
according to the energy considerations of their arrange- 
ments even at the edges of the systems. 

To obtain the pairing phase diagram we calculate the 
binding energy A as the difference in the ground state 
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energy Eqs when the on-site potentials are all off and 
after they were switched on: 

A(/(e)) = E GS (U A = U B = 0;/(£)) 

-E GS (u A ^U B ^0;f(e)), (3) 

where /(e) refers to the homogeneous, diatomic or Fi- 
bonacci arrangement of the atomic levels s. A positive 
A means there is local pair formation in the system. In 
Figures (2H3] we show the ground state phase diagram for 
binary Hubbard chains with homogeneous, diatomic and 
Fibonacci ordering of the site energies e. In all cases, 
the area above (and including) each line corresponds to 
paired states. The area below the line corresponds to 
non-paired states. 
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FIG. 4: Ground state pairing phase diagram for binary Hub- 
bard systems at half-filling. 
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FIG. 2: Ground state pairing phase diagram for binary Hub- 
bard systems with 2 electrons. FIG. 5: Minimum value of the attractive potential Ua for 

the systems to have electron pairing with a fixed repulsive 
potential Us/t = 8.0 as a function of the band-filling, n. 
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FIG. 3: Ground state pairing phase diagram for binary Hub- 
bard systems at quarter-filling. 

When browsing from Figure [D to Figure |H one ob- 
serves that the electron pairing is indeed enhanced due 
to quantum confinement and it also increases with the 



number of electrons in the system. In Figure [5] we com- 
pare our results directly to those obtained using Bethe 
Ansatz (Eq. [T]) , for the case of a binary Hubbard system 
with £a — £b — and two electrons. In Figure [21 it is 
also shown how a completely periodic modulation of e, as 
in the case of the diatomic chain, enhances the most the 
electron pairing, whereas an aperiodic modulation, such 
as that of the Fibonacci chain, enhances the electron pair- 
ing as well but this is only significant when compared to 
the homogeneous case. These results are comprehensible 
if we consider what was explained about the Fibonacci 
chains: There are more sites with energy values £,4 = 
(a total of 89) than sites with energy values Eb = 1 (a 
total of 55). This fact makes it plausible to observe the 
results of the Fibonacci chain rather closer to those of 
the homogeneous system than to those of the diatomic 
case. The same behavior of the phase diagram was found 
in the systems at quarter-filling (see Figure where it 
is even observed that the distance between the curves for 
both homogeneous and Fibonacci chains decreases. This 
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result indicates a sort of competition between these two 
systems as a function of the band-filling. An effect which 
close to half- filling renders interesting results. In Figure H] 
we show the results for the binary Hubbard chains at half- 
filling. While the diatomic chain remains well below the 
other two cases, the homogeneous and Fibonacci chains 
seem, at first glance, to behave the other way around, i.e., 
there are now more electrons paired for the homogeneous 
than for the Fibonacci system. We investigated further 
this behavior in the next way: After fixing a value of 
the repulsive on-site potential Ub, we obtained the min- 
imum value of attractive on-site potential Ua for which 
the systems have electron pairing. The results (see Fig- 
ure E|) show for the Fibonacci chain a uniform, though 
oscillating behavior in Ua- On the contrary, for the ho- 
mogeneous case, and for values of the band-filling close to 
n = 1, the minimum Ua values decrease faster than the 
Fibonacci case. This situation results in a critical value of 
the band-filling of n = 0.895 after which the behavior of 
the homogeneous and Fibonacci binary Hubbard chains 
is inverted. In general, quantum confinement enhances 
electron pairing in binary Hubbard chains. However, the 
topology of the confinement determines strongly the elec- 



tronic properties. The behavior of the homogeneous and 
Fibonacci chains close to half filling indicate that the 
quantum confinement effect is overridden by the effect of 
the aperiodic structure, as this is not happening in the 
diatomic chain. 

In conclusion, we considered a system of spin-i 
fermions in a one-dimensional lattice with position de- 
pendent atomic level £j, embedded in a binary on-site 
potential. Our results illustrate how the modulation of 
the atomic energy levels in a periodic and non periodic 
way affects the formation of local electron pairs in binary 
Hubbard chains with attractive and repulsive on-site po- 
tentials. We found a combined effect on the electron pair- 
ing due to the atomic level ordering and the increasing 
band-filling. The number of paired-states increases with 
the band-filling for the homogeneous, diatomic and Fi- 
bonacci chains. Moreover, there is a critical value of the 
band-filling for which the effect of aperiodicity compared 
to the homogeneous case is reversed. 
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